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Theoretical  studies  of  Fourier  telescopy  for  deep  space  imaging. 

Final  Report:  Contract  F61775-98-WE012 
Dr.  Peter  Bakut 

International  Informatization  Academy, 

Laser  and  Information  Technologies  Department 
Volokolamskoe  sh.,  95,  Moscow,  123424  Russia. 

September  1999 

1.  Introduction. 

In  paper  [1]  it  was  shown,  that  if  to  illuminate  a  target  by  a  pair  of  coherent 
emitters,  the  total  reflected  energy  is  proportional  with  the  Fourier  -  component  of  the 
target’s  image  on  spatial  frequency,  proportional  with  position  vector  between  emitters. 
Using  many  emitters  located  as  a  square  matrix,  and  measuring  a  total  reflected  energy 
from  each  pair,  it  is  possible  to  receive  a  good  many  the  Fourier  -  components  and  by 
inverse  Fourier  transform  to  restore  the  image  of  the  target.  In  the  same  paper  the 
method  of  phase  closure  was  offered  to  eliminate  unknowns  initial  phases  of  emitters 
and  to  save  the  problem  of  their  phasing.  This  method  of  imaging  of  remote  objects 
was  named  Fourier  -telescopy.  In  further  papers  [2-6]  this  method  was  more  deeply 
studied  and  was  perfected.  However  there  was  a  series  of  the  not  solved  problems  rather 
essential  to  understanding  of  the  Fourier  -  telescopy  method.  It  is  following  problems: 

As  the  total  reflected  energy  cannot  be  received,  what  happens  at  reception  of  a 
part  of  reflected  energy,  what  requirements  to  a  size  of  the  receiving  aperture? 

How  does  one  optimally  arrange  emitters,  that  it  was  possible  to  implement  a 
method  of  phase  closure  and  to  receive  the  greatest  number  the  Fourier  -  components 
at  minimum  number  of  emitters? 

How  do  a  partial  coherence,  frequency  drift  and  not  stability  of 
illumination  have  influence  on  Fourier-telescopy  process. 

What  quality  the  image  may  be  obtained  by  the  Fourier  -  telescopy,  what  is  its 
contrast,  speckle  structure,  resolution  of  details? 


0 1  -  54 7^ 
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These  problems  are  investigated  in  the  present  paper  and  the  some  answers  are 

obtained. 

2.  Statistical  model  of  a  field  in  the  Fourier  -  telescopy  including  micro 

roughness  of  object. 

In  the  present  section  there  is  calculated  the  field  on  the  receiving 
aperture.  It  is  taken  into  account  relief  and  micro  roughness  of  the  surface  of  the 
object. 

Let  emitters  are  located  in  points  p  q  ,  p  i  vj  P  •  Let's  designate: 

R  is  the  position  vector  of  some  central  point  of  the  target,  fi'om  which  we 
digitize  a  position  vector  of  surface  E ; 

n=  —  is  the  unit  vector  of  a  direction  on  target,  q  =  —n  is  the  unit  vector 
R 

of  the  direction  from  target  on  receiving  -  transmitting  position  (fig.  1).  The  field 
of  the  n—th  emitter  in  a  point  of  surface  of  the  target  is  equal 


u„  (r^ ,  m) « - — - g„  (vs  )u„  (m) , 

2%R 

where  gn(v^)=  (p)t/p  is  a  directivity  diagram  of  an  emitter, 

At 


(p)  is  the  field  in  aperture  of  an  emitter, 


RA- 
R  +  r-^ 


is  the  unit  vector  of 


a  direction  in  a  point  of  surface  of  the  target,  (fco)  is  a  spectrum  of  a 
signal,  emitted  by  n—th  emitter.  The  composite  field  of  all  emitters  is  equal 

- — - g„(vs)H„(to).  (2.1) 

„=o  2icR 

Further  we  suppose,  that  all  emitters  are  focussed  on  the  target  and  consequently 
the  directivity  diagrams  (v^  )  are  identical;  g'n  (v^  )  =  )  • 
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The  field  (2.1)  is  reflected  by  extended  object  with  surface  L,  on  which 
the  coefficient  of  reflectivity  is  distributed.  The  reflected  field  in  a  point 

p  of  the  receiving  aperture  is  equal: 

-{qn^  )K{rj^  )u^  (r^ ,  ,  (2.2) 


ik  1 


w,(r,/co)  =  — -  f- 
"  27t2.J 


R 


where  is  the  unit  vector  of  a  normal  line  to  surface  of  the  target  in  a  point  ^ , 
Substituting  (2.1)  in  (2.2),  we  obtain: 

-ik^-R-r-^  |-/A:|.R+r2  -p„  [ 


W. 


\2tiJ  2„^oj, 


X  iqnj^)K(r^)g(v^)u„(m)dL .  (2.3) 

The  obtained  formula  is  assume  for  further  reviewing. 

We  use  model  of  a  rough  surface  as  enough  smooth  underlying  surface,  on 
which  the  roughnesses  are  overlapped: 

r^=r  +  N^(r).  (2.4) 

where  ^  is  a  position  vector  of  rough  surface,  r  is  a  position  vector  of 

underlying  surface,  N  is  a  unit  vector  of  a  normal  line  to  underlying  surface, 
^(r)  is  height  of  roughness  above  underlying  surface  (fig.  2).  We  suppose, 
that  the  field  ^(r)  represents  homogeneous  on  a  underlying  surface  random  field 

with  zero  expectation  and  correlation  function  ^(^)^(^)  =  .S(^,^).  The 
homogeneity  means,  that  5(^5^)  =  5(/(fj,i^)),  where  /(^,^)  is  a  distance 

between  points  f[,  on  a  underlying  surface. 

Let’s  substitute  (2.4)  in  (2.3)  and  use  paraxial  approximation,  taking  into 
account,  that  R  »  p,  r,  p„ ,  i.e.  sizes  of  target,  matrixes  of  emitters  and  receiving 
aperture  are  small  in  comparison  with  distance  up  to  the  target.  Received  field  in 
point  p  of  receiving  aperture  can  be  noted  as: 
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where  =  qr,  =  ^p,  p„^  =  ^p„  are  projections  of  vectors 

f,  p,  p„  on  direction  q , 

f^=f-  q{qr)  is  perpendicular  to  q  component  of  vector  f , 

V„^(r)  =  q'V^(f),  V^(r)  =  p, - -  +  v - ,  u,v  are  local  orthogonal 

^  du  dv 

coordinates  on  a  underlying  surface  in  an  environ  to  a  point  F ,  p,  V  are  basis 

vectors  of  coordinate  axises. 

In  a  temporal  aspect  the  field  (2.5)  is  look  like: 

u,{p,t)  ~  2Re  =  2Ree'“o't/(p,0 , 

®0±Aco 

where  the  complex  amplitude  of  received  field  is  equal 


Here  is  complex  envelope  of  a  signal  illuminated  by  « —th  emitter.  As 

it  is  slow  function,  we  have  neglected  delays  on  rather  small  times 

.  The  delays  on 

2Rc  c  .  c 

if  the  target  is  on  a  vertical  (we  consider  the  emitters  are  located  in  a  horizontal 
plane).  However,  if  the  target  is  not  on  a  vertical,  there  are  additional  increments  of 
optical  path  from  different  emitters  up  to  the  target  and  from  the  target  up  to 
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different  points  of  the  receiving  aperture  at  the  inclination  of  a  sighting  line.  Let's 
while  suppose  signals  of  emitters  high-stable  to  neglect  delays  complex  envelope 


on  times 


Consequently,  we  have: 


-1  2ikR^  .  .X  r  +Mpg+Pnq+rg-Ngm) 


U,(p,t) 


2X^R^ 


K(r)g 

NqdS, 


(2.7) 


,  R 

Where  t  =t - . 

c 


3.  Reception  and  processing  of  a  reflected  field. 

In  the  present  section  the  reception  and  processing  of  a  field  Uj.{p,t) 
(2.7),  ended  by  deriving  of  input  data  for  elimination  of  unknowns  of  phases  by  the 
phase  closure  method,  deriving  the  Fourier  -  components  and  forming  the  image  is 

considered.  We  presiune,  that  the  amplitudes  (2-7)  have  personal 

frequency  coloring:  U„{t) all  differences  Aco„— are 
various,  also  bandwidth  of  signals  U^{t^  A©  « |Ao)„  —  Aco^j .  The  field 
U^(p,t)  transits  an  optical  filter,  is  detected,  after  detected  signal  transits  a 
system  of  filters  tuned  on  frequencies  A©„  —  A©^,  and  after  phase  detector  we 
obtain  primary  signals  (t) . 

In  [10]  the  problem  of  choice  of  frequencies  A©^  surveyed,  that  all 
differences  A©„  -  A©^ ,  n,m  =  0...N,  were  various.  For  basis  takes  minimum 


discrete  Aco  and  the  frequencies  are  searched  as:  Aco^=c(m)AQ.  For 
N  =  30,  for  example,  the  following  outcomes  are  obtained: 

Table  1. 


For  approximated  estimations  it  is  possible  to  note,  that  c(m)  ~  m  . 

In  outcome  of  photo  detecting  the  receiver  measures  power  of  a  field  (2.7) 
passed  over  the  receiving  aperture  A : 


^  m,n=\ 


2  _  2 

l(AC0„  —ik  —  'r2ik{p„q 

X(3.1 


)xjl 


2  2 
-ik^-^ 

e  ^ 


^11  2  I  ( Pm^ll  P« 2  ) 


,(p  7  ^^'^iA(r<il-N,xWx)y(r,2-Nq2W2)i , 

■yn~^L2)^ 


ss 

(  \ 

X  1-- 


q2  J 


where 


m  irv - 

«(4)  =  ^  dp 


is  a  point  spread  fimction  of  the  aperture  A . 

Passing  the  signal  (3.1)  through  a  system  of  filters,  tuned  on  frequencies 

Aco^  -  A(Of^  and  realizing  phase  detection,  we  can  separately  select  complex 
amplitudes 


2  _  2 

_il^Pjn_PjL+2ik{p.„  -p™ ) 

x^„(t)=u^{nu:{t')e 


F 

^  mn  » 


Pmn=  \  \e  *  *  «(>il-?±2)x 


where 
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X 


V 


N. 


J 


1 


V,^(r2) 


V 


N. 


ql  J 


xK{r,)K\r^)g 


(r  \ 

^11 

U  J 


^2ik\rqX  -Nq^  ^(Fj  ))-(rq2  -Nq2  ^(F2  ))] , 


*1  ^  \dr^^dr^2- 


R 


(3.4) 


Here  is  used  N^dS  =  df^,  and  S_i  is  silhouette  of  the  target  in  a  plane 
(perpendicular  target  sighting  line). 

The  signals  X^y^{t)  are  result  of  primary  processing  of  a  field.  The 


further  transformations  are  bound  with  digitization  of  signals  and 

computer  processing.  A  first  step  is  the  procedure  of  phase  closure  and  elimination 

2  _  2 

of  factors  {t')e  ,  containing  unknowns  phases.  In 

outcome  the  magnitudes  (3.4)  will  be  chosen  in  the  pure  state.  The  first  idea 

of  the  Fourier  -  telescopy  encompassed  that  the  magnitudes  are  the  Fourier  - 
components  of  the  image  of  the  target.  In  particular,  it  takes  place  at  the  rather 
large  receiving  aperture.  However  from  the  formula  (3.4)  is  not  seen  immediately 
Fourier  transforms  and  image  of  the  target.  Therefore  principal  problem  is  the  study 
of  requirements,  when  (3.4)  transmutes  in  a  Fourier  transform  of  some  image  and 
study  of  properties  of  this  image. 


4.  Case  of  the  rather  large  receiving  aperture. 

The  idea  of  the  Fourier  -  telescopy  method  consists  that  at  the  very  large 
receiving  aperture  A 


Thus  is  gave: 


F„„=X^R^  je 

Si 


...  (Pm-P«)4 


\  K  J 


5(4i-42) 


R 


N. 


\K{rt 


9  y 


g 


\Rj 


dr^  = 


(4.1) 
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=  T^R^F* 


r^(P/n  Pn) 


R 


(4.2) 


i.e.  there  is  a  Fourier  transform  of  function 


A(4)  = 


N 


i^(^)r 


9  ; 


\s 


\Rj 


(4.3) 


on  frequency  The  function  Ii(r±)  is  some  optical  image  of 

R 


the  target.  This  image  depends  on  a  reflection  factor  |Ar(r)|^,  intensity 


illumination 


g 


R 


,  and  from  geometry  of  a  surface  of  the  target 


N. 


Due  to  dependence  from  micro  roughness  ^(r)  this  image  is 


<1  ) 


random,  that  should  be  exhibited  in  it  spekleness.  However  this  dependence  is  faint, 
as  it  is  usual  a  steepness  of  micro  roughnesses  are  small  enough,  and 

V.^(r)  V„^(r) 

^  « 1  .  Neglecting  — - - ,  we  shall  receive: 


q 


N. 


/,(4)=|X(r)f 


\g 


f  f^\ 

\Rj 


(4.4) 


This  image  does  not  depend  at  all  on  geometry  of  a  surface  of  the  target. 

Let’s  clarify  now  as  far  as  be  great  the  aperture  A  owes,  that 

u{r^Y  —^12)  possible  to  consider  5—  figurative  in  relation  to  function 

under  the  sign  of  integral  in  (3.4).  Let's  mark,  that  length  of  the  spike  of 

function  “^2)  order  where  yfA—  diameter  of  the 

receiving  aperture..  Most  quickly  changed  factor  under  an  integral  (3.4)  is 
^ik^qx  ^q2  ^924(^2))]  function  did  not  vary  on  an  length 
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,  it  is  necessary  to  have  »  k^R,  ^JA  »  k^R  ,  where  kg,k^  are 

characteristic  steepness  of  a  underlying  surface  and,  accordingly,  of  micro 
roughness.  At  characteristic  steepness  of  roughness  about  0.1  [7]  it  reduces  in  silly 

large  values  of  diameter  of  the  receiving  aperture  (at  i?  ~  10  km  'JA  » 1  km'. 


XR 
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At  ~  40000  km  (geosynchronous  objects)  ^JA  »  4000  Atw).  It  is  put  under 
doubt  practical  reliability  of  the  Fourier  -  telescopy.  However  essential  role  in 
estimation  can  play  providing  for  influence  of  micro  roughness  of  object’s 
surface. 

5.  Consideration  of  deriving  of  distortionless  Fourier  coefficients  of  the  micro 

rough  target  image. 

Generally,  as  was  already  marked,  the  quantities  are  not  the  Fourier 

coefficients  of  any  target  image  ,  and  have  a  complex  form  (3.4).  As  into  the  formula 
(3.4)  go  into  a  micro  unevenness  i^(r),  quantities  are  random.  Therefore  it  is 

necessary  to  study  statistical  properties  of  quantities  and  to  clarify,  when  they 

become  look  like  the  Fourier  coefficients  of  a  image  of  the  target,  and  which  image. 

Let’s  calculate  statistical  characteristics  of  quantities  (3.4),  suggesting 

^(r)  by  a  homogeneous  field  on  underlying  surface.  Let's  mark,  that  for  such  field 
^(r)  and  V^(r)  independent.  Therefore 


7mn)=  j  ^  ^  «fel-42)x 


xKin)K\r2)g 


R 


g 


*  ^ 

I  R  J 


drndri_2- 


(5.1) 
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/  -2iA:[iV„i  ^(r,  )+iV„2  )1\ 

In  an  Apendix  it  is  shown,  that  (6  ^  )  represents  rather  narrow 


peak,  length  of  which  there  is  an  order 


Xh 


.  It  is  possible  to  expect,  that  this  peak 


8  —  figurative  in  reference  to  all  functions  under  the  sign  of  integral  in  (5.1).  The 

requirement  of  it  8  —  figurativeness  in  reference  to  function  C  ^  looks 

like: 

2n,  Lf«4TZ^R^R, 

R  4%g^ 

that  with  a  major  reserve  is  executed.  Here  is  a  maximum  size  of  a  matrix  of 

2ik{r  ’-’T  ) 

emitters.  The  requirement  of  8  —  figurativeness  in  relation  to  function  C 

gives  faint  limitation  in  a  reference  steepness  of  a  underlying  surface:  kg  « 'In  — — . 

k 

Remaining  functions  under  the  sign  of  integral  in  (5.1)  are  essential  more  slowly 
varied.  Thus,  we  have: 

{Fmn)  =  ^  \e 


R 


\m)f 


Si 


g 


u 


x<^ 


Al  %  > 


1^1 


A-^  \e 


V 


N, 


q2  ) 


)dri.2 


drn 


Zkg, 


\ 


.XP^-Pn)fL  ^1^2^ 


ik 


R 


n]  2al 


V  k  ) 


-*xi2 


\K(r)\ 


N.. 


g 


\r) 


dr^  «. 


Pn)^l 


\e  R  Io{r±)dr_i. 
%nGr  ^ 


(5.2) 


^  -^1 
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We  have  obtained  a  Fourier  transform  of  function 


Nl 


4(4)  =  ^ 


n]  2<yl 


C5\N 


1 


li  K  j 


^xi2 


(5.3) 


As  shown  in  [7],  /o(4.)  ^  optical  image  of  micro  rough  object,  which  is  obtained 

in  a  natural  light  with  the  help  of  a  standard  optical  instrument.  Let's  underline,  that  it 
differs  from  an  image  /i(^)  (5.4),  obtained  by  a  Fourier-telescope  at  the  super-large 
receiving  aperture. 

Let’s  investigate  now  fluctuation  of  .  For  this  purpose  we  calculate 


.  In  an  Apendix  the  relevant  evaluations  are  done  and  the  following 


outcome  is  obtained; 


2 

J  - 42)4(4i)4(42)^nL^^2i 

SlSi 


(5.4) 


We  should  be  interested  with  relative  quantity  of  fluctuation  of  : 


^mn 


Using  (5.2)  and  (5.4),  we  obtain: 


11 


mn 


J  J^(4i  ”  42  )4  (4i  )4  (42  y4i^42 

_ 


\e  ^  4(4¥4 


2 


(5.5) 
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For  an  estimation  of  this  quantity  we  shall  assume  at  first,  that  the  aperture  A  is  small 
also  sizes  of  the  target  less  of  a  resolution  element  « - .  In  this  case 

*  A 

h{r^l  —  f^2)  ^ 'Hmn  matters  about  unites..  It  evidences  about  strong 
fluctuations  ,  so  the  concrete  realization  ofF^„is  random  also  unpredictable.  Let's 
assume,  that  the  aperture  A  is  large  and  there  are  many  elements  of  a  resolution  are 
stacked  on  the  target: 

o  2  2 

M  «  '  - . » 1 ,  i.e.  A  »  ^  ^  .  (5.6) 

Thus  ^(^1  -  r_i2  )^R  A,  A’ S  -  rj_2  )  and  we  obtain: 


^mn 


Sl 


!■ 


R  /o(rj_)JrL 


(5.7) 


7^  R^ 

From  here  it  is  visible,  that  at  A  » -  r|^„  « 1 .  It  means,  that  practically  do 

Sfg 


not  fluctuate  and  are  equaled  medial:  Thus,  at  execution  of  a 

requirement  (5.6)  output  signals  of  a  Fourier-telescope  represent  Fourier- 


components  of  an  image  of  the  target,  this  image  /o(Pl)  expressed  by  the  formula 
(5.3)  and  represents  the  usual  image  obtained  with  the  help  of  a  standard  optical 
instrument  in  a  natural  light. 


6.  Reception  of  a  field  with  matrix  of  spaced  telescopes. 
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It  was  above  established,  that  the  size  of  the  receiving  aperture  of  a  telescope  for 
forming  images  should  obey  to  usual  Rayleigh  measure: 

— —  »  Stg  ,OT  M  =  y  »  1  , 

where  M  is  number  of  a  resolution  elements,  stacked  in  square  of  the  target.  However 
for  usual  telescopes  the  requirement  of  high  optical  quality  of  a  receiving  mirror  adds  to 
this  requirement.  In  the  Fourier  -  telescopy  the  receiving  telescope  appears  as  a  collector 
of  an  energy,  and  any  key  requirements  to  quality  of  a  mirror  is  not  showed.  Therefore 
for  the  Fourier  -  telescopy  rather  actual  is  represented  idea  to  synthesize  the  large 
receiving  aperture  from  the  small  spaced  apertures.  The  problem  of  a  precise  coherent 
bringing  together  of  light  beams  from  the  spaced  telescopes  here  does  not  exist. 

Let’s  assume,  that  the  illuminated  power  receives  by  K  telescopes  with  the 

identical  entering  apertures  Aq,  placed  on  terrain  in  points  i9„,  n  =  1,...,/^ .  Let's 

designate  a  pupil  functions  of  telescopes  A{p  —  •  Let's  note  a  received  field  (2.7) 

in  the  reduced  form; 
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,  where  Mq  —  “  I  ^  ^  ^  point  spread  function  of  the  aperture  , 

y  y 

A) 


S.  ik 


r^k 


G(f)=  is  a  lattice  factor  of  receiving  telescopes.  If  receiving  telescopes 

it=l 

spaced  in  nodes  of  a  square  lattice  with  a  pitch  dj. ,  then 


■wy' 

A  - rUh 

A:=l 


(6.2) 


A/? 


where  —  an  integer  vector.  The  function  (6.2)  is  periodic  with  a  period 

df. 

Lj.  in  two  directions.  At  enough  major  K  it  looks  like  a  narrow  peak,  concentrated 


IZ  .  ( kif\\  -  f\y ) 

above  a  site  cr^  =  — ^ .  The  function  Wq  - 

K 


V 


|g(^1  —  r_i2  )  is  equivalent  to 

f2 


a  point  spread  function  of  the  uniform  major  aperture  if 


A 


2  and  L^.  exceeds 


enough  size  of  the  target  (that  imder  the  sign  of  integral  (6.1)  there  was  only  one  peak 
of  function  Gi^^A  “^2)  )•  Thus,  if  the  linear  size  of  the  target  is  having  filled 

AJR. 

square  A  by  small  telescopes  with  distance  between  them  dj,  < - we  obtain  for  the 

^tg 

Fourier  telescopy  the  same  effect,  as  an  energy  is  received  by  one  major  telescope 
with  the  aperture  A . 

7.  Method  of  phase  closure  and  influence  of  temporal  unstability  of  illumination. 


The  method  of  phase  closure  allows  to  eliminate  of  unknowns  initial  phases  of 
emitters  and  thus  to  save  a  problem  of  their  special  phasing.  For  implementation  of  a 
method  of  phase  closure  it  is  necessary,  that  the  adjacent  emitters  were  shifted  from  each 


15 


Other  on  d  or  on  an  axes  jc ,  or  on  an  axes  y  (not  simultaneously).  The  examples  of 
configurations  of  matrixes  of  emitters  are  shown  in  a  fig.  4. 

Let  a  signal  (3.3)  from  ffz-  th  and  n-th  emitters  is  selected.  Let's 

consider  that  m<n  and  that  path  from  emitter  to  emitter  p„  pass  through 
emitters  P;„+i,P;„+2  vjPn-1  •  Simultaneously  signals  from  neighboring  emitters 
^{k-\)ki^)’  A:  =  m  +  are  received.  As  it  was  shown  in  Point  3  signals 

have  form: 


qn  yqm  > 


mn  ’ 


where  are  expressed  by  the  formulas  (3.4)  and  under  particular  requirements 

represent  Fourier  transforms  of  the  image  of  the  target.  The  phase  closure  method 


consists  in  a  multiplication  of  all 


^{k-\)k  (0 
^{k-\)k  (0 


A:  =  m  + 1,...,  n,  and  closing  of 


the  obtaining  product  by  multiplication  on  {t) .  Thus  is  obtained 


”  X 

r„„(0=  n 


{k-\)k  (0  _ 

mn 


k=m-{-l  \^(k-l)k  (0 


•  ^  2/r  ^  . 

k=m^l  ^ 


=  Uminunn 

After  temporary  smoothing  we  obtain: 

n  X 


f  In  ^  X} 

^1 


^mn  ~  ~ 


^idk-dk_{) 


F 


In 


\ 


{cipi  ctj^ ) 

j 


,  (7.1) 


where 


®(/)  =  argF(7),  F„„ 


‘V n,(t'W‘n(t-) 

w 


dt. 


(7.2) 
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As  a.]^  —  a]^_i  =  (±1,0)  or  (0,±1)  only  two  phases 


a  =  0 


f  2% 


L  271  ^ 


^,0  ,  (3  =  0  0,^ 


V  ^  J 


are  presence  in  expression  (7.1).  Let's  assume 


V  ^  y 


y  =  (a,  (3) .  It  is  easy  to  see,  that 


k=m+l 

From  obtained  expressions  follows; 


=m+l  \  ^  J 


(7.3) 


Y  =V  e 

^mn  ^ mny 


-iy(a„-an) 


F\  y(4-««) 


(7.4) 


If  to  put  =  V  =  const ,  the  quantities  represent  in  the  pure  state  Fourier- 
components  with  a  phase,  proportional  frequency.  At  image  reconstruction  we  receive: 


A  p 


2%p 


V 


2%i  _ 

—p{rL  -Y) 


(7.5) 


Thus,  the  restored  image  appears  shifted  on  a  vector  y,  that  is  unessential  for  the 
research  of  an  inner  pattern  problems  and  of  the  image  recognition. 

It  is  essential  to  mark,  that  as  a  result  of  phase  closure  phase  fluctuations  of 
signals  were  purely  eliminated.  It  is  valid  for  high-stable  signals,  space  length  of 
which  is  well  exceed  the  size  of  target.  In  particular  it  relate  to  chirp  with 
parameters  1 0  MHz  /  (X  sec . 

As  it  is  visible  from  expression  (7.4),  amplitude  fluctuations  of  a  signal  at  phase 
closure  completely  are  not  cancelled,  that  is  exhibited  available  of  random  factors 


at  the  Fourier  coefficients.  However  in  a  case  XjA/’  »1,  where  is  a  illumination 
pulse  duration,  the  temporal  smoothing  is  equivalent  to  a  statistical  average  and  we 
obtain;  ^  ^  const ,  i.e.  the  amplitude  stability  does  not  influence  on  the 

Fourier  -  teleseopy  process. 

Generally  in  case  of  partial  coherent  signals  the  output  signals  X^^^it), 
obtained  by  Fourier-telescope,  will  look  like: 
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2  _  2 


Fmnit)’ 


M  /•  -;/c - 

Where  J  Je  ^ 


_  -Pnhz) 


^  w(4i-rx2) 


X 


A 


X 


1 


N. 


c 


) 


z  \\ 


1 


N. 


q2  ) 


^2ik\rq^  -Ng^  ))-(rg2  -Ng2  Wl ))]  ^ 


XV 


m 


f'  _  ^(pg  P«g  ^gl)l  ,* \^,  _  ^(pg  P/wg  "h  ^^2) 


x^(ri)r(r2k[^^ 


A  Jr^^\ 


g 


c 

''12 


V  ^  y 


dr±idr±2  • 


(7.6) 


Here  functions  U^{t')  represent  a  slow  part  of  amplitude  and 


,,  ^{Pq+Pnq+^qlf 


—  quick  fluctuations  .  Usually  v(/)  approximate  by 


CT 


stationary  stochastic  process  with  time  of  correlation  T .  Thus  l^=  —  there  is  a  length 

2 


of  a  coherence  of  illumination. 


The  functions  v„ 


..  2(3 


q  Pnq  ^gl> 


and  v» 


Pq  Pmq  ^q2 , 


cannot  be  taken  out  from  under  the  sign  of  an  integral  because  of  presence  of  a  delay 


2r„ 


Most  important  in  the  practical  ration  there  is  the  case,  when  A/qX  » 1,  where 


c 

A/’o  is  a  bandwidth  of  filters  selecting  frequency  C0^„ .  Thus  the  smoothing  realized 
after  phase  detection,  practically  will  reduce  to  average  of  a  signal.  If  the  signals,  emitted 


different  emitters  are  independent,  we  receive:  {ti  )v„  (^2  )  ®  I^'Vnfe)=0. 

From  here  follows,  that  between  signals  emitted  by  the  different  emitters,  rigid 
correlation  should  exist.  It  is  possible  in  case  the  emitters  feed  from  one  generator  and 
the  optical  paths  from  the  generator  to  emitters  are  sustained  identical  within  the 
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accuracy  of  lengths  of  a  coherence  of  the  generator.  Thus  we  have: 
Vm  (0  **  (0  =  "'^(0  ^Iso  we  obtain: 


-  -  ^  —IK - t-lK -  .  V 

Pn.„(0-‘f„„(0=  j  je  '>  ^ 


5,  -S' 


±^1 


X 


rr^^ 


1  ^ 


vAcn) 


N. 


J 


1  ^ 


V,^(r2) 


N. 


q2  ) 


,2jl[(r^l  -Nqx  ))-(rq2  -Ng2  ^(^2  ))] 


X 


XV 


m 


~ P«g  ^(p^  Pmq  ^ql) 

t  I 


xK(r,)K\r2)g 


/r  \ 
^±1 


yRy 


g\^\dr±idr^2  = 


R 


=  I  Je  « 

S± 


^  w(4i-42)x 


r 


X 


1  ^ 


N„ 


\  ; 
Xlj/f 


1  ^ 


V,^(r2) 


N. 


^liklrq,  -Nq,  \{r,  ))-(rq2  -Nq2  m  ))]  , 


Pw^r  Pmq  ^q\  ^q2 


q2  ) 

\ 


V 


K{r^)K\r^)g 


J 


KR) 


S\  dr±idA2- 
.  K  y 


(7.6) 


For  physical  estimations  we  approximate  correlation  function  by 


gaussian  function:  l|/(A^)»^  .Thus we  obtain: 


FmniO^e 


2 

ip mg  ~Png  ) 
21} 


II- 

iSj_  Sj_ 


R 


X 


V 


N. 


rr 


J  V,^(r2) 


V 


N. 


2.^-Pmg  P;,g 


^  w(4i-rL2)' 


2 

(>gl->V2)-^""*~y--Kl^(n)+^g2^(f2)) 


2/, 


X 


q2  ) 
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xK{r,)K\r^)g 


fr^,\ 


\Rj 


g 


^12 


R 


dr^^dr^2 


(7.7) 


Now  we  conduct  an  average  on  ^(r)  according  to  a  procedure  explained  in  an 
Appendix.  Let's  obtain: 

iPmq-Pnq)^ 


mn 


it)) 


in 


R 


STCGt 


i 


1  + 


liNi 


{2kl,f  a|  N] 


\+i 


.Pmq  Pnq 


“i2 


2k 


_ 

1  iM'H^q  ^ 


2  _2»r2 


2 


^.i2 


l^(^)| 


N. 


g 


UJ 


dr. 


If  the  target  is  on  a  vertical  axes,  so  =  0 ,  we  obtain: 


(7.8) 


r 


iN 

■■A-  ^ 


Sira 


1' 


ik 

e  R 


^iVj_ 


1  + 


2  9  ^ 

i  Kj 


1 


\^{2kl,) 


2^ 

llNl 


X 


2{ki,y 


2  ^£^9 

1+(2W,)2^^ 

A/". 

.xe  ^  ^ 


|^(^)|' 


N„ 


g 


/■* 

rj_ 

Cr) 


dr. 


From  (7.9)  it  is  visible,  that  the  image  will  register  in  this  case 


(7.9) 


2{kn) 


i'o(fj=e 


■  'M 


k^l 


n  2  9  ^ 

1} 


9  J 


i 


CT?V? 


i^(?)r 


1  UeiV^ 

l  +  (2iL)2^ 


g 


R 


llNl 


(7.10) 


20 


This  image  differs  from  /qC^)  particular,  it  depends  on  length  of  a  coherence 

of  radiation  4  •  At  major  4  it  go  to  Iq  (r^) .  And,  4  should  not  be  especially  major: 

^4  » ~ ~ lOj  l^»X. 

k 

However  it  is  very  strong  on  a  form  of  the  obtained  image  renders  inclination  of  a 
sight  line  in  a  case,  when  the  target  is  not  on  a  vertical.  As  it  is  visible  from  expression 
(7.8)  the  signal  from  a  pair  of  emitters,  the  difference  of  optical  paths  from  which  in  a 
direction  to  the  target  is  comparable  to  length  of  a  coherence  of  radiation.,  fades. 
Therefore  angle  of  deflection  of  an  optical  axes  from  a  vertical  should  obey  to  an 

inequality:  Sina  «  — ,  where  Lf  -  a  maximum  size  of  a  lattice  of  emitters. 


8.  Image  reconstruction  and  estimation  of  its  quality. 


The  secondary  signal  processing  in  the  Fourier  -  telescopy  consists  in  image 
reconstruction  on  signals  .  Let  requirements  formulated  in  the  previous  section  are 

carried  out,  and  are  proportional  with  the  Fourier  -  components  of  image  /(4.) 
of  the  target: 


mn 


\lirje 


.j.iPm-Pn)ri 
R 


dr=F 


(  KPm-Pn)^ 

R  J 


(8.1) 


Let's  suggest,  that  the  emitters  are  disposed  in  nodes  of  a  rectangular  lattice  with  step  d , 
i.e.  in  points  with  position  vectors  ,  where  a^,  m  =  0,. . . ,  iV,  are  integer  vectors. 
Then  the  quantities  proportional  with  the  Fourier  -  components  on  frequencies 

2nd 

multiple  to - ,  where 

L 

X-/  — 

d 


(8.2) 
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The  area  of  an  integration  in  (8.1)  natural  restricted  to  sizes  of  the  target.  Formally  it  is 
possible  to  take  it  of  any  greater  size.  If  distance  between  emitters  d  is  small  enough, 
thus  a  maximum  linear  size  of  the  target  <  Z,  then  area  of  an  integration  in  (8.1) 

can  be  restricted  to  quadrate  Lx  L.  Then  it  is  possible  on  coefficients  (8.1)  with  the 
help  of  a  Fourier  series  to  restore  the  image  periodically  repeating  on  a  plane  with  period 
L  in  two  directions  without  superposition  of  periodically  repeating  duplicates..  The 
accuracy  of  restoring  is  defined  by  number  used  the  Fourier  -components.  Let  set  of 


(A/^  +  1)^  quantities  F 


^jPm  P«) 

R  J 

/  2np^ 

L  J 


contains  Fourier-components 

,  p  gG  , 


(8.3) 


where  G  is  some  two-dimensional  set  of  integer  points.  The  restored  image  looks  like: 


^  peG  k 


(8.4) 


J 


The  degree  of  proximity  of  functions  /  (F)  and  /(F)  characterizes  a  resolution 
capability  of  a  Fourier-telescope  defining  Fourier-components  (8.3).  It  is  easy  to 
see,  that 

7(F)  =  JjG(F  -  r')I{r')dr' ,  (8.5) 


where 


2ni _ 

L  p^Q 


The  function  G(r)  looks  like  a  peak,  concentrated  above  a  site  by  a  size  a  =  — ,  where 

M 


M  -  number  of  points  in  a  set  G ,  i.e.  number  the  Fourier  -  components  used  for  image 
reconstruction.  From  here  follows,  that  the  linear  resolution  capability  of  Fourier- 
telescope  6  is  equal: 


5« 


L 

4m' 


(8.7) 
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Let’s  remark  now,  that  the  set  of  {N  + 1)  of  quantities 


F  -F 

^  mn  ^ 


^(Pw  Pn) 
R 


contains  different  Fourier-components  F 


,  and  many 


from  them  are  doubled  including  multiply.  As  any  process  of  measuring  is  accompanied 


by  noise,  the  output  signals  of  Fourier-telescope  ~  F 


K9m-9n) 

R 


are  measured 


with  errors.  From  here  follows,  that  the  backup  of  the  Fourier  coefficients  should  be  used 
for  increasing  of  accuracy  of  image  reconstruction.  Taking  into  account  expression  (5.2) 


we  can  write: 


Snai  * 

^mn  ~  j2^2  A  ~  ^ 

LX  A 


KPm  -P») 
R 


^mn  > 


where  measurement  errors  we  consider  independent  and  Gaussian.  From  this  the 
following  algorithm  for  optimum  image  reconstruction  on  output  signals  can 


be  obtained: 


m,n=0 


(Pm-PnFl 


^  dr,  =  min . 


Varying  on  I{r±),  we  obtain  the  equation: 


fr/-  \ 


m,n-0 


-(“w-S/i)?! 


where  5^  is  integer  position  vector  of  w  —  th  emitter.  It  is  possible  to  enter 


characteristic  function  of  m  —  th  emitter 


^ot]’ 


representing  infinite  matrix  with 


unity  in  a  point  of  standing  of  m  -  th  emitter  and  zero  in  remaining  points. 
Equation  (8.8)  will  be  copied  as: 


j7(4,)EA(*)e’ 


dfn=  .(8.9) 

W2,«  =  0 
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where  A(^)  =  —  5^]  (8.10) 

m=0 

is  a  characteristic  function  of  all  transmitting  matrix  of  emitters  representing  an  infinite 
matrix  with  unites  in  points  of  standing  of  emitters  and  zero  in  remaining  points. 
Conversing  (8.9)  by  the  Fourier,  we  receive: 

2ni  __  ^ 

where 

C{p)  =  ^KhKk-p)=  Y.^^-^mW-P-'^n]=  YAPm-^n-pY^-^^) 

is  a  number  of  pairs  of  emitters,  the  vector  of  distance  between  which  is  equal  p .  The 

I 

sense  of  expression  (8.11)  is  very  prime:  the  total  of  quantities  ,  representing  the 

,  is  divided  on  total  nxunber  of  such  ,  i.e.  for 

image  reconstruction  the  average  Fourier  coefficients  are  used. 

From  (8.1 1)  the  following  algorithm  of  image  reconstruction  implies: 

1  —pf,  1  N 

KA)  =  7^  Z  ^  ^  77^  YKrA^m  -^n-pl  (8-13) 

P  psG  ^\P)m,n=0 

Let's  calculate  statistical  properties  of  the  restored  image.  First  calculate  average  value. 

Using  (5.2)  it  is  easy  to  receive: 

27tf  jy 

(^(4))  =  72  Z  fcfe-««-p]=4(4)-(8-i4) 

P  peG  ^\P)m,n=0 

Thus,  the  recovered  image  on  the  average  is  equal  /qC^i)  (3-3). 

Let’s  calculate  correlation  function  of  the  restored  image.  We  have: 


same  Fourier  coefficient  F 
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2ni, 


~  '  x(8.15) 

^  p,q^G 


1 


C(p)C(q)^,„%^o 


1  ^iPhi-qhi)  1 

-  Z  — 


t  [{Knni)-{Kn){Hl)hm-an-pW^ 


JU  n 


p,qeG 


C{p)C{q)^^^^^^A 


J  J  KKl  ~  ^12  Vo  (^ll  )^0  (^12  )  ^ 


^  -^b^-5i)rii-(5„-a,)r[2Lr 


m,n,kJ-0 


^[Gm-Gn-pWk-Gl-q] 


\drndrl2 


J  i4g(4i  -;ii)G'(42  -?;2)|Wu  -^Lf 


5,  5 


1  *^1 


N  _^5  F 

Ujnr 


where  H(f)=  ^ 

m=0 


We  shall  consider  a  case,  when  the  aperture  Aj.  is  great  enough  and  it  is  possible 
to  consider  h{ArjJ  as  5—  figurative  in  relation  to  functions  under  the  integral  (8.15). 

It  is  easy  to  note,  it  is  necessary  that  diameter  of  the  receiving  aperture 

considerably  exceed  a  maximum  linear  size  of  a  matrix  of  emitters  :  -sfA  »  Lf .  In 

this  case: 


.2 


where 


^  P 


(8.16) 


(8.17) 


From  here  we  discover  a  relative  variance  of  fluctuation  of  luminosity  of  the  restored 
image  (contrast): 


g(2i,4) .  _  Li 

aP  ~A 


(8.18) 
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Lr 

At  —  « 1  C  « 1 ,  i.e.  the  image  not  fluctuates  and  is  equal  to  average  value. 

A 

Except  for  contrast  of  the  important  property  of  the  image  is  the  medial  size  of 
spots  .  The  medial  square  of  a  spot  can  be  evaluated  under  the  formula: 


^ _ 


WnAx)  L] 


4^’ 


where  6  = -  is  element  of  a  resolution. 

A 


(8.19) 


In  paper  [7]  the  uniform  measure  of  quality  of  the  image  in  a  coherent  light 
equals  to  a  relative  variance  of  correlation  of  the  medial  image  and  its  random 
realization  was  entered: 


({j  l{HAWA))B{rnA2  F  ) 

{j(^(4))'^4F 


(8.20) 


The  substitution  in  this  formula  of  expressions  (8.13)  for  (/(^j_))  and  (8.14)  for 
^AlAl)  gives: 


AS,^ 


(8.21) 


Where  M 


is  number  of  elements  of  a  resolution  on  the  target. 


Let’s  consider  the  problem  on  an  optimum  configuration  of  emitters.  From  a 
disposition  of  emitters  it  is  required,  that  at  their  given  number  the  greatest  number  the 
Fourier  -  components  should  be  obtained.  For  matching  different  configurations  by  this 
measure  we  use  the  formula  (8.12),  which  gives  number  the  Fourier  -  components  with 

firequency  p : 
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C{p)  =  '^KhK^-p)-  (8-20) 

ic 

The  formula  (8.20)  shows,  that  number  the  Fourier  -  components  with  frequency  p 
equal  to  number  of  self-intersections  of  a  configuration  of  emitters  at  its  shifted  on  a 
vector  p .  In  fig.  4  the  variants  of  filling  of  a  fi’equency  plane  for  configurations 
figured  in  fig.  5  are  figured.  In  the  Table  2  the  amounts  of  emitters  for  each  of  < 
surveyed  configurations  providing  the  same  number  M  of  the  Fourier  coefficients  are 
represents. 


Table  2. 


Configuration 

amount  of  emitters 

Y-figurative 

~L3Vm 

L-figurative 

~L4Vm 

Il-figurative 

^  1.5a/ Af 

0  -  figurative 

~1.8Vm 

+  Figurative 

~2Vm 

From  the  table  it  is  visible,  that  the  Y-configuration  is  most  favorable,  and 
the  -i--configuration  is  least  favorable.  For  providing  of  the  same  quality  of  image 
reconstruction  it  is  required  in  Y-configurations  1.5  as  many  emitters  as  in 
+  -configuration. 


9.  Inference  and  conclusions. 

In  paper  the  following  outcomes  are  obtained. 

1.  The  quality  of  the  image  obtained  by  the  Fourier  -  telescopy  method  is  determined  by 
two  factors:  by  a  size  of  the  receiving  aperture  both  size  and  number  of  emitters  of  a 
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transmitting  matrix.  For  deriving  the  images  with  given  number  of  elements  of  a 
resolution  the  size  of  the  receiving  aperture  A  should  obey  to  usual  Rayleigh  measure: 


»  Sfg 


I.e.  the  necessary  number  of  Rayleigh  elements  of  a  resolution  should  be  stacked  on 
object.  However  in  a  case  the  Fourier  -  telescopy  the  receiving  telescope  appears  as  a 
collector  of  an  energy  and  there  is  not  showed  any  special  requirements  to  optical 
quality  of  a  mirror  of  a  telescope  .  In  particular  a  large  receiving  aperture  can  be 
implemented  with  the  help  of  spaced  small  telescopes.  The  small  telescopes  should 


place  in  nodes  of  a  grid  with  a  pitch  dj.  < - and  fill  m  square  A . 

2.  The  transmitting  matrix  can  be  implemented  as  two  or  several  bars  of  emitters 


located  in  nodes  of  a  grid  with  a  pitch  d  < - .  For  a  possibility  of  embodying  of  a 

^tg 

phase  closure  method  each  folio  wing  emitter  of  bar  should  be  disposed  in  adjacent  on  x 
or  y  a  node  of  a  grid.  The  element  of  a  resolution  at  image  reconstruction  equal  to 


T^R^  2 

5  «  — - — ,  where  M  is  number  of  different  Fourier-components.  M  ~  N  ,  where 

d^M 

N  —  number  of  emitters.  The  configuration  of  a  matrix  of  emitters  should  be  such  to 
provide  maximum  number  of  different  Fourier-components.  Were  explored  Y-,  L-,  I-,  O- 
and  +  figurative  configiuration.  The  greatest  number  of  Fourier-coefficients  has  appeared 


for  Y-figurative  configuration 


,1.3, 


and  least  for  +-  figurative 


3.  At  embodying  the  formulated  requirements  on  sizes  of  receiving  and 
transmitting  of  the  apertures  it  is  possible  to  receive  the  image  given  by  expression  (5.3): 
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Such  image  receives  by  a  standard  optical  instrument  in  a  natural  light. 
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4.  The  temporal  stability  of  illumination  practically  does  not  influence  on  quality 
of  the  recovered  image.  The  slow  temporal  stability  is  cancelled  as  a  result  of  the 
procedure  of  phase  closure,  the  fast  temporal  stability,  bound  with  a  partial  coherence  of 
illumination,  averages  at  signal  processing. 

In  this  paper  we  have  not  touched  problems,  bound  with  moving  of  the  target 
and  presence  of  Doppler  shift  of  frequency.  Besides  the  work  of  a  Fourier-telescope 
on  the  remote  space  target,  as  a  rule,  carries  at  low  levels  of  signals  in  conditins  of  the 
score  of  quantum.  A  research  of  process  the  Fourier  -  telescopy  in  these  requirements 
represents  major  interest.  These  problems  nor  have  found  reflectings  in  this  report.  It  is 
desirable  to  prolong  work  on  the  Fourier  -  telescopy  in  these  directions.  At  last,  the 
theoretical  researches  are  bound  to  making  of  mathematical  models  frequently  requiring 
of  experimental  inspection.  Making  an  experimental  bench  on  the  Fourier  -  telescopy  and 
holding  of  laboratory  experiments  therefore  is  desirable. 
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Apendix. 

Calculation  of  ^ .  Agrees  (4.1): 
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Since  ^(F)  is  Gaussian  we  have; 

-  2ik{N^,m  -  N2,m))\  ^  g- " 2Nl,N2,b(ri,r,)) 


^-4kV^Nl^(\-bir[,r2)) 


(A.2) 


B(r  Y  ) 

Here  is  designated:  af  =  B(f, r),  b(fi , ^2 )  = - ^  .  We  have  taken  advantage  of 

B{r,r) 

that  radius  of  correlation  of  micro  uneveimesses  ^(r )  is  rather  small  in  comparison  with 


reference  scales  of  a  relief  of  a  underlying  surface.  Let's  suggest  that  »  A-.  In  this 

assumption  the  function  (A.2)  is  rather  acute  the  peak  concentrated  in  limits  of  locally 
flat  site  of  a  underlying  surface.  In  limits  of  this  site  it  is  possible  to  enter  orthogonal 
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coordinate  system  m,v  .  We  use  for  analytical  estimations  Gaussian  approximation  of 
correlation  function: 


(til  -Vif  +{U2-V2f 

21} 


b{^{uy,v^),r{u2,V2))^e 

Thus  the  function  (A.2)  can  be  approximated  as  follows: 


-  2ik[N^^{r{u^,v{))- N2qUf{u2,V2)%,  ^  p-  -  e 
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Supposing  r  =  rQ  +  p,*M  +  V-  V,we  can  obtain: 

((a^i,  -  V,«r,))K  -  + 

a^(n) , , ..  9|(a) 

9«i  ■'  5v,  3«2  «  5v2 
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(A.5) 


Let's  mark,  that  the  function  (A.4)  is  6 -figurative  in  relation  to  function  (A.4). 
Character  length  of  a  peak  of  function  (A.5),  agrees  (A.3),  makes  .  Character  length  of 
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/c  O'  e 

a  peak  of  function  (A.4)  makes  - - ,  that  at  — ^  » 1  is  much  less.  We  suggest 

X. 

also,  that  mikro  unevenness  are  not  resolved  by  the  aperture  A ,  i.e.  the  function  (A.3)  is 
8  —  figurative  in  relation  to  function 

Using  made  estimations,  we  obtain: 

-V,4(P2)))rfS2  » 
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Substituting  (A. 6)  in  (A.l)  we  obtain: 

StIOi: 


(A.6) 


(A.7) 


Calculation  ) .  We  have: 
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r  \  fr 
^±1  *  'L2 

R)  KR 

Having  taken  advantage  that  )  is  Gaussian,  we  obtain: 

/g-  2ik[N^fi(r{)  -  +  N,^^(r,)) 


^m)K\r^)K\^)K(r,)g 
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g  dS^dS2dS^dS4 
\  R  J 


^  NlqJ^ZqKh-h)-  Vl,JV3,6('i.'3)  +  ■V„A?4,6(?i.?4))  ^ 

^  g-4i'{V3,V3,6(?2,?3)-  Vj,V4,6(f2.?4)  -  N,^N q.Kh,^))  g, 

This  function  has  two  peaks  with  single  maxims  on  diversities  ^  ^  ^  and 

^  ^ ^  •  These  peaks  are  rather  acute,  practically  5  —  figurative.  Therefore 

function  (A.8)  can  be  represented  as  the  total  of  these  peaks.  For  anal3dical  estimations  it 
can  be  used  Gaussian  approximation  of  peaks  similar  (A.  3) 

For  approximating  of  a  peak  in  an  environ  of  a  diversity  ^  ^ ^  =  F4 , 

let's  assume  r2  =  Fj  +  3c,  F4  =  +  5^ .  Let's  suggest,  that  the  peak  concentrated  in  the 

field  of  small  x  and  y  in  limits  locally  of  flat  sites  of  a  underlying  surface  in  an  environ 


of  points  ^ ^ .  Let's  receive: 

/(Fi,Fi  +F,F3,F3  +;;)  = 

-Ak^(5\  Nlq  (1  -  Z)(Fi ,  Fi  +  F))  +  Nlq  (l  -  bir^, ,  F3  +  y)) 

—  Q.  X 

^  Ak^^lN^N^q  [-  b{r^ ,  F3  )  +  b{r^ ,  F3  +  >))  +  b{r^  +  F,  F3  )  -  b{f^  +  x,  F3  +  y)] 
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(A.  10) 


For  approximating  of  a  peak  in  an  environ  of  a  diversity  ^  ^  j  ^2  “  ^4 » 

assume  r4=^+35.  The  peak  concentrated  in  the  field  of  small 

X  and  y  in  limits  of  locally  flat  sites  of  a  underlying  surface  in  an  environ  of 
points  ^ ^ .  Let's  obtain: 
f(7i,72’7i  +x,72  +y)  = 
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The  substitution  (A.12)  in  (A. 8)  gives: 
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We  here  have  taken  advantage  of  the  formula: 
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Using  the  formula  (5.3),  we  can  write: 
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Following  outcome  here  is  used: 
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It  is  easy  to  note,  that  second  term  in  (A.  14)  neglectfully  small  in  comparison  with  the 
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first  term.  Their  ration  varies  from - ;r  at  the  small  apertures  A  up  to  — .  .  .  at  the 


large  apertures.  Therefore  finally  we  have: 
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In  particular  we  obtain: 
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PD-photodetector, 

Z-  summator, 

CO  -bandpass  filter  tuned  to  friquency  C0„^, 
LPF,  LPFi  -  low  pass  filtrs. 

Fig.  3 


42 


Fig.5b.  Fillness  of  friquency  plane 
at  LI -configuration 


Fig.  5b.  Fillness  of  friquency  plane 
at  L-configuration 


Fig.Sd.  Fillness  of  friquency  plane  at 
0-configuration 


